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I. INTRODUCTION 


As long as a full quantum theory of gravity is not available, it is important to address the 
quantization of particular models. Among the most interesting ones are systems of spheri¬ 
cally symmetric gravity plus matter, because they are sufficiently non-trivial and nevertheless 
simple enough for some technical treatment. Moreover, they fold physical applications in 
the study of gravitational collapse and quantum aspects of black holes. 

In this paper we discuss the Lemaitre-Tolman-Bondi (LTB) model which describes a 
self-gravitating dust cloud. It has originally been introduced by Lemaitre [lj in order to 
study cosmology where it has indeed found interesting applications, cf. [2] for details and 
references. Here the focus is on the canonical (Hamiltonian) formalism for both the classical 
and the quantum LTB model. The classical part is mainly intended as a preparation for the 
quantum model but it exhibits interesting aspects onits own. 

The physical questions that one eventually seeks to address include the issue of singularity 
avoidance in the quantum theory, the quantum evolution of black holes, and the role of the 
naked classical singularities in quantum gravity. In the simpler case of a single dust shell , 
for example, it was shown by an explicit construction that a unitary quantum theory exists 
in which the classical black-hole singularity is avoided, cf. SB- In fact, a wave packet 
describing a classical shell evolves into a superposition of collapsing and expanding shell and 
yields destructive interference at the place of the classical singularity. The formalism of loop 
quantum gravity has also been applied to spherically symmetric systems and it has been 
claimed that the singularities are avoided 0. 

In recent years, earlier works on the canonical quantization of the Schwarzschild black hole 
0,0 have been applied to develop a canonical description of the collapse of a marginally 
bound, spherical, time-like dust cloud 0 0. ‘Marginally bound’ means that the various 
shells of the dust cloud start with vanishing velocity at infinity. Quantization of this classical 
dust system leads to the Wheeler-DeWitt equation for the wave functional describing the 
quantum collapse. As special cases, it was possible to obtain the Bekenstein mass spent rum 
and statistical entropy of the charged and uncharged black hole 


10]. Later, in 


11 


it 

was shown that the semiclassical (WKB) treatment of the Schwarzschild black hole in this 


canonical picture describes Hawking’s thermal radiation and in 12], going beyond the WKB 


approximation, that quantum corrections render the spectrum of the radiation non-thermal. 

Here we again apply quantum geometrodynamics and the Wheeler-DeWitt equation to 
the LTB model in order to see how far this approach can be developed. In contrast to 
the earlier works, we consider the generic case, that is, the case including the non-marginal 
models for which the classical shells start with a non-vanishing velocity at infinity. This paper 
is organized as follows: in Sec. II we introduce the LTB model and present the canonical 
formalism for the general models. Some technical details are relegated to Appendix A. In 
Sec. Ill we perform a lattice regularization and fold for a particular set of factor orderings 
exact solutions to both the Wheeler-DeWitt equation and the diffeomorphism constraint. 
The uniqueness of this set of solutions for the given ansatz is shown in Appendix B. Sec. IV 
presents a brief summary and an outlook on possible future work. 
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II. THE CLASSICAL LTB MODEL 


A. Metric and classical solutions 


The LTB model describes a self-gravitating dust cloud. Its energy-momentum tensor 
reads T /w = e(r, p)u fl u u , where = w M (r, p) is the four-velocity vector of a dust particle 
with proper time r and labeled by p ( p thus labels the various shells that together form the 
dust cloud). The line element for the LTB spacetime is given by 

ds 2 = -dr 2 + ^ 9pR \ . dp 2 + R\p){d6 2 + sin 2 0d0 2 ) . (1) 

1 + 2 E(p) 

Inserting this expression into the Einstein equations leads to 

8nGe(r, p) = R and (M) 2 = ~r + 2E ' ( 2 ) 

where F(p) is some non-negative function, G is the gravitational constant, which we set 
equal to one in the following (except for some key expressions where it is retained), and we 
set c = 1 throughout. The case of collapse is described by d T R(r, p) < 0. 

There exists still the freedom to rescale the shell index p. This can be fixed by demanding 

R(Q,p)=P, ( 3 ) 


so that for r = 0 the label coordinate p is equal to the curvature radius R. Now we can 
express the functions F(p) and E(p) in terms of the energy density e at r = 0. From (12]1 
one gets 

F(p) = 8tt f e(0, p) p 2 dp , (4) 

Jo 

1 Ayr CP 

E(p) = ~[d T R(T = 0,p)} 2 -/ e(0, p)p 2 dp . (5) 

z P Jo 

The interpretation of these quantities is that F(p)/2 is the active gravitating mass inside 
of R(r,p), while E(p ) is the total energy of the shell labeled by p. The marginally bound 
models are defined by E(p) = 0. In the present paper we discuss the general case which 
includes the non-marginal case defined by E(p) ^ 0. 

The solution of (J2J) is given by 


T ~ T 0 (P) = 


RHt,p)Q (- HMlp) 


v'W) 


where r 0 is a constant of integration that can be fixed by (j3j) to read 

p 3 ^Q(- 2 -¥) 


to (p) = 


Vf 


( 6 ) 


(7) 
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The function Q(y) is defined by the expressions 


Q(y) = | 

( arcsin (y/y) 


{ y 2 

Q(y) = 3 for y = o , 

Q(y) = | 

I arcsinh (yf—y) 

l (—2/)* 


V 1 -y \ 

V ) 


for 1 > y > 0 , 


v 71 -y \ 

y J 


for 0 > y > —oo , 


( 8 ) 

(9) 

( 10 ) 


where Q > 0 and 1 > y > —oo. 

Equation (J6J) shows that at the dust proper time r = ro(p) the shell labeled by p has 
reached a curvature radius R = 0, that is, the physical singularity. So r can only take values 
between —oo and r 0 (p). 


B. Hamiltonian formalism 


For the canonical formalism one starts with the general ansatz for a spherically-symmetric 
line element, 


ds 2 = -N 2 dt 2 + L 2 (dr + N r dt) 2 + R 2 dQ 2 , 


( 11 ) 


where N and N r are the lapse and shift function, respectively. The canonical momenta are 
given by 


P L = § + NrR ‘ 

P R = — \-LR -LR+ ( N r LR )' 


N L 


( 12 ) 

(13) 


A dot denotes a derivative with respect to coordinate time t, while a prime denotes a 
derivative with respect to r. All variables are functions of t and r. 

A Legendre transformation from the Einstein-Hilbert action leads to 


Seh = J dt J dr ( P l L + P r R - NH 9 - N r H 9 ) + S d x , 


where the Hamiltonian and the diffeomorphism (momentum) constraint are given by 

PlPr 


H 9 = -G 


R 


Hf) + i 

2 R 2 G 


L R' 2 / RR' \'' 

~2~2l + 


H 9 = R'P r - LP' l , 


(14) 

(15) 

(16) 


respectively, and the boundary action S'as is discussed below. 

The total action is the sum of (TT4|) and an action S d describing the dust. The canonical 
formalism for the latter was developed in 13], cf. also [ii| . It reads 


S d = 


dt J dr (P t t - NH d - N r H d ) , 


(17) 
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where the Hamiltonian and momentum constraints are 


H a = P T \ 1 + 


r 


/2 


L 2 


and Hi = r'/'. 


( 18 ) 


In principle one would prefer to take a fundamental held (e.g. a scalar held) for the matter 
part. However, this would make the formalism much less tractable [14]. Moreover, the main 
features discussed here already exhibit themselves for the dust model. 


C. Mass function in terms of canonical variables 


In the following we shall write the mass function F(p ), which was introduced in (|4[), 
from the canonical data. This is essential for deriving consistent falloff conditions that are 
appropriate for a realistic collapse model. For the Schwarzschild case this was done in [7], 
while for the marginal solutions (E = 0) this was done in [9]. Here we shall extend the 
derivation to the non-marginal case. 

We start by demanding that the spacetime described by the metric (TTT]) be embedded in 
a LTB spacetime. Considering the LTB metric (JT[) , a foliation described by functions r(r, f) 
and p(r, t ) leads to 


ds 2 = 


r 


p\d P Rf 
1 + 2 E 

p'\S p Rf 


if _ 1 2f _ mMi , didr 


1 + 2 E 


(19) 


t ' 2 ] dr" + R 2 dfi . 


1 + 2 E 

Comparison with (TTTT) gives the following three equations: 

p\d p Rf 
1 + 2 E 


N 


L 2 N r2 = f 2 


L 2 N r = 


Pp\d p Rf 
1 + 2 E 


- 2 _ p'\d P RY 


L = 


— TT 


J 2 


Eliminating N r from 


1 + 2 E 

with the help of (JUD and (J22D gives 
/V 2 = (0„Rf {fyf - tV)2 


L 2 (l + 2E) ’ 

Dehning E =1 — F/R we can rewrite Einstein’s equation 

(■ d r R ) 2 = 2E + 1 - T . 


as 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 


Now we insert the expressions for lapse function (1231) and shift vector (12T1) into the expression 
for the canonical momentum P L (fT2l) . Dehning R = (d p R) /\/l + 2 E and J = p'f — pr' we 
have 


-R- ~ R7 + R'[-fr' + pp'R~ 


(25) 
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Using the expression (1221) for L on the right-hand side as well as 

R = fd r R + pdpR = fd T R + RpV 1 + 2 E , (26) 

R! = t' 8 t R + p'd p R = t' 8 t R + RpW 1 + 2 E , (27) 


we get 


LP l 

R 


RJ = -R 2 Jp d T R - RV 1 + 2 Et'J , 


and therefore 


LP l ( R' - t'8 t R) 8 t R 


R 


VT+2E 


y/1 + 2Et' 


0 R'\/2E + 1 — F Ft’ 


Solving this equation for t' gives 

1 


t' — —— ( R'Vl + 2E -P + 


VI + 2 E Vl + 2 E 


LP l VT+2E 


R 


Inserting this expression into (1221) yields 


t 2 /2 /2 

L = R p — r 


R 


/2 


1 + 2 E 1 + 2 E 


2 R’t’ , ---- Pr' 2 

Vl + 2 E-T - 


1 + 2 E 


= J IR 


L 2 P 2 


?/2 


which leads to 


R^_ Pi 

^ ~ ~l 2 Ip ' 


We can thus express F locally in terms of the canonical data as follows: 


(28) 

(29) 

(30) 

(31) 

(32) 

(33) 

(34) 

(35) 


F 


R 


Pi 

1-)-h _ 

R 2 


R' 2 

U 


(36) 


This is the same expression as was obtained in |9] for the marginal models. It thus possesses 
a much wider range of applicability and holds, in fact, for all cases. 

Since R = F at the horizon, F = 0 there. We can check that though T appears in the 
denominator of (1311) , t 1 is well behaved at the horizon, as it should be: 


T ’ i=» i (g + L) . (37) 

As in the case of the Schwarzschild black hole [j] and the marginal LTB-model [q|, one 
can make a canonical transformation in order to elevate the mass function F to a canon¬ 
ical coordinate. The expressions are the same as in these earlier papers. The canonical 
transformation is 


(r, i?, L, P T , P R , P L 


(r, P, F, P T , P R , P F ) 


(38) 
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where 


( 39 ) 



A = (RR')(LP l )' - (RR!)'(LP l ) . 


(40) 


The action in the new canonical variables then reads 


/ 


f 

Jo 


oo 


d t 


dr P T t + P r R + P f F -NH- N r H r ) + S&e , 


(41) 


where the new constraints are 



+ 4 GPP f P r ) +P t a/1 + —, 


(42) 

(43) 


//, = t ' p t + + / '/v. 


We shall now discuss the boundary action Sqt, in more detail. 

D. Boundary action 

Boundary terms are obtained from a careful discussion of the falloff conditions for the 
canonical variables. This was investigated for the marginal case in |!f. For the non-marginal 
case, this derivation remains unchanged if E —» 0 for r —► 0. But this follows if one assumes 
that the initial density profile is regular at the center: Taking R( 0, r) = r (this is just a 
choice of scaling), then 



(44) 


where /i(r) = e(0,r) is some initial density profile, cf. ((4j). For a profile that is regular at 
the center, choose 


OO 



n =0 


which gives 


OO 



n=0 


As the center of the collapsing cloud is taken to be at rest in a spherically symmetric collapse, 


the velocity profile d T R(0,r), must, to leading order, behave as some positive power of the 


label coordinate, r. It then follows from the ‘velocity equation’ (1241) that E behaves as r 2 


or higher. In particular, E —> 0 as r —> 0. 

The only boundary term is obtained from the variation with respect to L and reads 



(47) 


where N + (t) = N(t,r oo) is the lapse function at infinity and M + (t ) = F(r oo)/2 is 
the ADM mass. To avoid the conclusion that N + (t) is constraint to vanish, which would 
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freeze the evolution at infinity, the boundary term has to be canceled by an appropriate 
boundary action. This can be achieved by adding the surface action 

Sas = - J cIt N + (t)M + (t) . (48) 

Since varying N + would lead to zero mass, Kuchar has argued in that N + has to be treated 
as a prescribed function. The lapse function gives the ratio of proper time to coordinate 
time in the direction normal to the foliation. Since N r (r ) vanishes for r —>■ oo, the time 
evolution at infinity is generated along the world lines of observers with r = const. If we 
introduce the proper time, r + , of these observers as a new variable, we can express the lapse 
function in the form N + (t ) = f+(t). This leads to 


S aE = -'JdtM + f + . (49) 

Thus we have removed the necessity of fixing the lapse function at infinity. (In [3] this is 
called ‘parametrization at infinities’.) 

In |S|, the proper time f+ was identified with the dust proper time at infinity, t+. By 
realizing that f+ is a priori not directly related to t+, we propose here a different treatment. 
Extending the treatment from the Schwarzschild case [7] to here, the aim is to cast the 
homogeneous part of the action into Liouvilee form and to find a transformation to new 
canonical variables that absorb the boundary terms. This can be done by introducing 
the mass density T = F' as a new canonical variable and using that F( 0) = 0 (which is 
appropriate for a collapse situation). Part of the Liouvillc form can then be rewritten as 
follows: 


6= / d rP F 8F - M.8f, 


drP F (r) / d r'ST(r l )+f + 6M + -S{M+f + ) 


= / d rP F {r) / d r , 8T(r , ) + '^- / d r'5Y{r') - 5{M + f+) 

Jo Jo 2 J o 

roc ro o roc — 

= / drhT(r) / d r'P F (r')+ / d r'5T(r')—- 5(M+f + ) 

Jo Jr Jo 2 

/ o° / — roc \ 

dr <5T (■+ J d r'P F (r') \ — S(M + t + ) . 

In the calculation we have used the identity 

roc rr roc roc 

/ d rP F (r) / dr' hT(r') = / drhT(r) / d r'P F {r') , 

Jo Jo Jo Jr 

which can be seen by integrating 

pr rr \ f rr roo 

dr'/V'V) / dr'dr(r') ) = P F (r) / dr'ff(r') - S\ tr) / dr'/Vlr'j 


(50) 

(51) 

(52) 

(53) 

(54) 

(55) 

(56) 
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from 0 to oo. The left-hand side gives zero and the rest the desired equation. 

From line (1541) we can read off that Pr = f+/ 2 + J r °° drPp. Thus we see that Pr(oo) = 
t + /2. As will become clear from the next subsection, this means that the Killing time at 
infinity matches the prescribed function f + . Thus the new action reads 


S'eh = j dt J d r[P T f + P r R + P r T - Nil" - N r H" r 


The constraints in the new variables are 


H = — 


1 fYFf 


//,. = t'P ; 


2 L \GF 
R'Pr 


-AGFP'Pr )+P t \I 1 + 


r f - 


L 2 


yp;. . 


(57) 

(58) 

(59) 


The Hamiltonian constraint can be greatly simplified if the momentum constraint is used 
to eliminate Pr = —Pr. The details of the calculation can be found in Appendix A. The 
constraints (1581) and (1551) can then be replaced by the following equivalent set, 


H = G ( P T 2 + PP 2 r ) - 


4 GT 


0, 


H r = r P T + R'P r - TP r ' « 0 . 


(60) 

(61) 


These equations will be used as the starting point for the quantization in Sec. III. 

We emphasize that the relative sign between the dust kinetic term and the gravitational 
kinetic term can change because T > 0 (< 0) outside (inside) the horizon. Since this change 
of sign is already present in (1581) . it was not introduced by using the momentum constraint 
to eliminate Pr. In fact, this change of sign of the gravitational part is also present in the 
corresponding expression for the Schwarzschild black hole, cf. Eq. (117) in [7j. Thus, it 
arises due to the choice of the new canonical variables, which are obtained by embedding the 
spatial spherically symmetric hypersurface into a LTB or Schwarzschild spacetime. To define 
F, we had to use the canonical momentum Pr. Thus we have mixed the original canonical 
coordinates with their momenta, which means that the new configuration space, spanned 
by r, R , and F, is different from the original one containing the three-geometries, which is 
spanned by r, R, and L. The change of sign is of fundamental interest in the quantum theory 
because the original Wheeler-DeWitt (WDW) equation has a (locally) hyperbolic kinetic 
term, which is of importance for the formulation of the proper boundary value problem 
3, 15]. A change of sign has hitherto been noticed for the WDW equation in the presence 
of a non-minimally coupled scalar held 16 


E. Relation between dust proper time and Killing time 


In this subsection we explain under which circumstances it makes sense to interpret Pr 
as Killing time. Equation (j5Tl) is given in the new variables by the expression 

Rl 

t' = 2P/V1 + 2 E + —Vl + 2 E-F . (62) 
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Defining a = l/\/l + 2 E gives 


2 P' V 1-EE 

a ciE 


( 63 ) 


If the mass density vanishes for all r greater than a given r b , and if E is constant for all r 
greater than r b , Eq. (1631) can be integrated. This yields 


ar = 2Pr + / d R 


Vl - a 2 E 
E 


= 2P r + F\ D—+ ln 
1 - a 2 /2 


1 - Vl - a 2 E 


1 + Vl- a 2 E 


Vl — 


In 


a z 


\/l — a 2 P — Vl — a 2 


\/l — a 2 P + a /1 — 


(64) 

(65) 


The result can also be written in the form 

'Vl-a 2 P 


ar = 2Pr + F 


1-E 


2 tanh 1 Vl — a 2 P + 


Vl - 


: tanh 


-i vl 


Vl - a 2 E. 


( 66 ) 


Here we have assumed that 0 < a < 1 (corresponding to E > 0). Eq. (1241) guarantees that 
the discriminant 1 — a 2 E is non-negative. This condition may be violated in the quantum 
theory, so it is of interest to give the result for all cases. For a > 1 but R < Fa 2 /(a 2 — 1), 
the result can be found by analytic continuation to read 


ar = 2Pr + F 


Vl - a 2 E 
l-E 


+ In 


1 - Vl - a?E 


1 - Vl - a 2 E 


2 — a 2 \J a? — 1 

H- . tan 


a/ a 2 


Vl - a?E_ 


(67) 


Another analytic continuation then gives the result in the region R > Fa 2 /(a 2 — 1), 


a 2 — 1 
a 2 - 1 

( 68 ) 

These expressions will be of relevance in the quantum theory, see Sec. III. 

We know from Birkhoff’s theorem that the spacetime around a collapsing dust cloud is 
given by the Schwarzschild solution. In [l] it was shown for the Schwarzschild geometry 
that 2Pr is equal to the Killing time T. Thus the equation derived above connects the dust 
proper time with the Killing time at the boundary r b . For small T and E' the relationship 
may still be used, since then we have a small amount of dust propagating in a Schwarzschild 
background; only in this case the concept of Killing time makes sense. In the limit a —> 1 
(P —> 0) we obtain 


ar = 2P r + iF 


V a?E - 1 „ , —- (1 — a 2 /2) 

^---2 tan” 1 V a 2 E - 1 + v , ! In 

1 -E V^l 


V d 2 E - 1 - V 


V o?E - 1 + V 


r = T + 2VP 


^-A ln (y±y) 

2 l rR-yr/ 


(69) 


which is identical to the relation used in 


11 


for the marginal case. The plus sign after the 


T on the right-hand side has been chosen in order to describe a collapsing dust cloud. 
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As has been discussed for example in |17j, Eq. (T64|) gives for a Schwarzschild spacetime a 
relation between Killing time and the time used by families of freely falling observers. Each 
family is characterized by a fixed value of E. All observers within one family start at infinity 
with the speed Voo, where E = r;^ 0 /(2(l — v^)). In the marginal case they thus start with 
zero initial speed. In this case Eq. (IPOD gives the relation between the Schwarzschild time 
and the Painleve-Gullstrand time [171 ] . 


F. Reconstructing E and tq from the canonical data 


The aim of this subsection is to reconstruct the local energy E and the singularity curve 
To from the canonical data. We have performed this reconstruction already for the mass 
distribution function, see (136|) . The three functions E, r 0 , and F determine the classical 
LTB model completely. There are several reasons for reconstructing these quantities. First 
of all, they have a clear physical meaning, whereas this is not so evident for the canonical 
coordinates. This is the reason why we have replaced the canonical coordinates L and Pp 
with F and Pp. One might hope that it is possible to turn E or r 0 into canonical coordinates, 
although this has not been achieved yet. But even without this, they help interpreting the 
canonical coordinates. 

We shall first demonstrate the following connection between the energy densities T and 


Pr- 


Pr 


T 

2V2E + 1 


(70) 


for an arbitrary foliation. Consider the momentum constraint 


/ P'Pk 

p F r 


(71) 

°- T+ 

■ + 

Pr 


/ & / 

r 2 

P F T 

(72) 

= '-W 

4 P? T + 

P ’ 

1 T 


where, in the second step, we have used the Hamiltonian constraint (l60j) in the form 


Pr = - 


Pr 

T 



(73) 


Comparing (1721) with (1621) we can read off Eq. (1701) . which is also a constraint. Equation 
(1701) then gives us a simple expression for E in terms of the canonical variables, 


1 + 2 E 


r 2 m r 2 
4P 2 " 4 EPl - r 2 /(E) ' 


(74) 


Note that, using the solution ([6]) of Einstein’s equation, we have an expression for To — r, 
that is, for the remaining proper dust time until the dust shell reaches the singularity, 


r 0 - r = 


R Wq (—2|S) 

VF 


(75) 
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Since we know now how to express E by canonical coordinates, this equation gives us an 
expression for the bang time To in terms of canonical coordinates. 

Now one can express Pr in terms of F and E. Inserting the expression for P T into (1731) . 
we have 


T 2V2E + 1 

1 r 

F 2V2E + i 


(±1)V2E + 1 - T 
<), R . 


(76) 

(77) 


G. Hamiltonian equations of motion 


Here we shall give the Hamilton equations of motion and derive Einstein’s equation 
from them. The Hamiltonian equations are generally given by 


X = {X,n[N]+7i r [N r ]}, 

Px = {Px,H[N} + H r [N r ]} , 

where we have introduced the smeared constraints 

/»oo /»oo 

H[N}= / d rN(r)H(r), H r [N r ] = / d rN r (r)H r (r) 

Jo Jo 

Starting from the action (1571) . the Hamiltonian equations of motion are 1 
r = 2NP t + N r r' , 

P T = ( N r p r y , 

R = 2 NFPr + N r R' , 


Pr = ~N 


FP, 


T 2 F \ 
R 2 ' AF 2 R 2 ) 


R 


+ 


+ (N r P R y , 


i = (ATT)' , 


P r = N — + N r PJ + 
1 2F r 


d fN(f) ( + 


T 2 (r) 


R(f) 4JF 2 (r)R(f) 


Consider now 


R N r =o N2FP r jtoJ TPr2\/2E + 1 

= -N2P7 = -r-’ 

which we can solve for Pr- Inserting this expression into the Hamiltonian constraint 
and again using (1701) gives 


(78) 

(79) 


(80) 

(81) 

(82) 

(83) 

(84) 

(85) 

( 86 ) 

(87) 


0 = H = 


r 2 


+ 


r 2 (d T Fty- 


r 2 


4(2 E + 1) T4,(2E +1) 4T ’ 


( 88 ) 


1 Note that 6F(r)/5T(r) = 9(r — r). 
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Solving for ( d r R ) 2 leads to Einstein’s equation (J2]). Note that we did not have to specify 
the lapse function. 

One can show, moreover, that F, E, and tq are constants of motion, that is, they have 
vanishing Poisson brackets with the Hamiltonian constraint. In their given local form they 
do, of course, not commute with the momentum constraint, since the latter generates their 
transformation with respect to a relabelling r —* f(r). One would, however, expect that 
a suitable non-local form commutes with H r and thus turns them into real ‘observables’ 
(similar to what one would expect to happen with the geometric operators in loop quantum 
gravity 0]). This would coincide with the interpretation of them being the physically relevant 
variables energy and bang time. 

We also note that because of (j8Tj) . N has no longer (for vanishing shift) the interpretation 
of being the ratio of proper time to ADM time, cf. [4j. The reason is that we have squared 
the original version of the Hamiltonian constraint, see App. A. If we define a new version of 
the constraint by taking a square root, 

H : =Pr- « 0 . (89) 

we find that 


{r,W t [Af>]} = Af> (90) 

and recover for the lapse function N^ the old interpretation. 

We finally remark that the algebra of the constraints cannot be of the general form (given 
for example in |4|), because we have used the momentum constraint to eliminate Pp in the 
Hamiltonian constraint. In fact, a short calculation gives 

{H[N\,H[M)} = 0, (91) 

{H, r [AT] , H[N]} = H[N, r N r - NN{ , (92) 

{H r [N r l Hr[M r ]} = Hr [[AT, M r ]\ . (93) 

We note that the Poisson bracket of the Hamiltonian with itself vanishes, Eq. (1911) . in 
contrast to the general case where it closes on the momentum constraint. The other brackets 
coincide with the general case. The transformations generated by the Hamiltonian constraint 
can thus no longer be interpreted as hypersurface deformations. They are in general not 
orthogonal to the hypersurfaces, but act along the flow lines of dust. 


III. DIFFEOMORPHISM INVARIANT QUANTUM STATES 
A. Quantum constraints 

We shall now apply the quantization procedure proposed by Dirac and turn the classical 
constraints into quantum operators, cf. |4|. Starting point are thus the expressions (1601) and 

m- 
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The translation of Poisson brackets into commutators is achieved in the Schrodinger 
representation by substituting 


Pr{r) 


h 5 
i 8r(r) ’ 


P R (r) 


h 8 

i SR(r ) ’ 


Pr(r) 


h 5 

7 <sr(r) 


(94) 


and acting with them on wave functionals. The Hamiltonian constraint (1601) then leads to 
the WDW equation, 



+ F «|y + ^ F)i(0) 


5 

8R(r) 

$[r(r') I fi(r') i r(r')] = 0 1 


(95) 


where A and B are smooth functions of R and F that encapsulate the factor ordering 
ambiguities. We have introduced divergent quantities such as 5(0) in order to indicate that 
the factor ordering problem is unsolved and can be dealt with only after some suitable 


regularization has been performed, cf. 18]. That is, one would like to choose the terms 


proportional to 5(0) in such a way that the constraint algebra closes, which is usually called 
‘Dirac consistency’. 

Quantizing the momentum constraint (16TT) by using (1941) gives 


R'- 


5r(r) 5R(r) 


-r 


5T(r) 


T [r(r / ),i?(r , ),r(r')] = 0 . 


(96) 


Up to now, the quantum constraint equations have been formulated only in a formal way. 
The next subsection is devoted to the application of a lattice regularization. 


B. Lattice regularization 


We follow here the suggestion made in 12] and consider a one-dimensional lattice given by 
a discrete set of points r l separated by a distance a. In order that the momentum constraint 
is fulfilled in the continuum limit, it is important to start with a corresponding ansatz for 
the wave functional before putting it on the lattice. We therefore make the ansatz 


T [r(r), R{r), T(r)] = U[ dr T(r)W(r(r), R(r), T(r)) 


(97) 


where U : M —> C is at this stage some arbitrary (differentiable) function. Using T in the 
exponent instead of R' or t' is suggested by the form of the WDW equation (absence of 
derivatives with respect to T) and the fact that F' = T is related to the energy density. 
The ansatz has to be compatible with the lattice, which means that it has to factorize into 
different functions for each lattice point. So we have to make the choice U = exp, which 
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gives 


where 


T[r(r), R(r), T(r)] 

(98) 

= exp drT(r)W(r(r), R(r), F(r))\ 

(99) 

= exp ^lim ^ ^ r tWi (r(r*), i^fa), F(^))j 

(100) 

= lim TT exp (cr Tj Wj (rfa), i?(r;), F(r*))) 

cr—>0 

i 

(101) 

= lim TT T,; (r(r•*), Rfa), Tfa), F(r *)) , 

cr—>0 

i 

(102) 

i 

F(r ,) = £ <xr,. . 

(103) 


As 


m 


j =o 

12] we implement the formal expression 5(0) onto the lattice as follows, 


5(0) —> lim — 

o-^o a 


The lattice version of the WDW equation (t95|) then reads, 


Gti 


f d 2 d 2 d \ 

(^ + ^ + W 3 o^) 


+ B(R,F) + 


a 2 r 

4GT, 


= 0 • 


(104) 


(105) 


We now insert the ansatz (11011') and make for convenience the redefinition W = iW/ 2. This 
leads to 


cr 2 T 2 


2 (dW(r, R, F)' 


Gh 2 i — 1WA2 - L ] + G h 2 F ( dW ^ 


dr 

aT 


GF 




limit 


+ 5(i?,F) = 0. (106) 

In order for this to be fulfilled independent of the choice of a (and thus also in the 
a —> 0) one is led to the following three equations, 

- dW(r, R, m 
OR 


Gh 


dW (r, R, T) 
frr 

{ d 2 


+ F 


FF [Gh 
d 2 




. dr 2 <9i? 2 


+ V«, r) Jj)M-'( T , fl, r) = o, 


(107) 

(108) 


and 


B(R,T) = 0 . 


(109) 
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The first equation, (1 1071) . the Hamilton-Jacobi equation, is the same as in [9], Eq. (5.11) 
and in [ll|], Eq. (12). The second equation presents an additional restriction on solutions 
of (1 1071) . The last equation (11091) tells us that working on the lattice is only possible if the 
factor ordering does not contribute to the potential term. If we find solutions to all three 
equations, we can do all other calculations on the lattice, since these solutions have a well 
defined continuum limit and satisfy the momentum constraint. 

At this point we have to make a few comments on the regularization procedure we are 
using. It was already noted that the lattice regularization does not solve the factor ordering 
problem. The lattice regularization just represents an ad hoc regularization in which the 
divergent terms have to cancel each other. Put differently, it is equivalent to a DeWitt 
type of regularization (which means setting J(0) = 0) with an additional constraint on the 
solutions. It has to be noted that the Hamilton-Jacobi equation. Eq. (11071) . can equivalently 
be obtained from the highest order of a WKB approximation 


11 


We have already emphasized that the signature in the kinetic part of the Hamiltonian 
constraint (160|) can change from elliptic (outside the horizon) to hyperbolic (inside the hori¬ 
zon). This thus occurs for the kinetic term of the WDW equation (1951) . too. As alr eady 
noted above, a similar phenomenon can be found in the context of quantum cosmology 


16 


As discussed in 19], we can say that the part inside the horizon is always classically allowed, 
whereas this is not necessarily the case for the outside part. The usual initial value problem 
appropriate for hyperbolic equations can thus only be applied for the region corresponding 
to the black hole interior. 

Since we have found from our special ansatz the two equations (11071) and (11081) . we cannot 
expect for them a well posed initial value problem to hold because the system is in general 
over determined. The goal pursued in this paper is not the discussion of boundary value 
problems but to find a class of exact solutions to all quantum constraints and to draw from 
them physical conclusions. 

It is instructive to check that one can get Einstein’s equation ([2]) from the Hamilton- 
Jacobi equation (1 1071) . A class of solutions is given by the complete integral 

.VT 


Gh W (r, R, a, b) = b + ar ± / d R- 


a 2 F 


T 

which is identical to the expression occurring in (1641) . cf. also Eq. (3.3) in [1 
made an ansatz of the form T = exp (iS/Gh) = exp (i f (F'/ 2) W), the relation 

F' d 

P R = - W 

R 2 dR 

F’ y/1 - a 2 F 

~2 F 

should hold. Using the Hamiltonian equation of motion for zero shift we have 

R = 2NFP r . 

Inserting this in (11121) gives 



= 1 — a 2 + a 2 


F 

R 


( 110 ) 
Since we 

(mi 

( 112 ) 

(113) 

(114) 
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This equation should be equivalent to Einstein’s equation given by 


d R 
dr 


= f + 2£ - 


( 115 ) 


For a = 0 the right-hand side of (111411 is equal to 1. This is a limiting case. As discussed in 
fl7| , it corresponds to observers that start at using the Eddington-Finkelstein coordinate 
v as their time coordinate. For an arbitrary a^Owe can write 


1 F 
~ a 2 + R 


(116) 


and can thus identify 2E = 1/a 2 — 1. Hence, a has the same interpretation as in (1641) . 
Letting a go to zero implies that E diverges. 


C. Solutions 


We have obtained above the two equations (110711 and (I108|l . which have to be satisfied 
in order to get a diffeomorphism invariant solution to the WDW equation. (We set now 
G — 1 — h.) Looking for particular solutions of the separating form W = ck(t) + (3(R), 
we recognize immediately that this system of equations is consistent only for special factor 


orderings. The particular factor ordering A(R,F) = F / {2R 2 ) used in 
among them. (For the WKB solutions used in 11] this is irrelevant.) 


11 


and 


3 


is not 


Tackling the problem from the opposite point of view, one can ask for which factor 
orderings we do get a separating solution. We find 


A(R, F) — 


2 R 2 


1 + 


1 -a 2 T 


(117) 


This leads to 


W(t, r, R,a) 


const. ± ar ± 



Vl - a 2 F 
F 


(118) 


These are identical to the solutions that were obtained in the solution of the Hamilton- 
Jacobi equation, cf. (1 11 011 and (166T) . Thus we can again identify 2 E = 1/a 2 — 1. Since 
classically E > —1/2, it follows that a should be real. Using ([TOjl . we can show that this is 
consistent: 


F t 'b a = = 


2\/2E + 1 




where we have dehned 


* 0 [r,fl,r] = exp(- / dr FW(t. R, F.a) 


(119) 


( 120 ) 


The integral appearing in (II 18(1 has been evaluated above for the various cases, see f[66ll . 
([67]1 . (l68|l . We recognize in particular from (16811 that the wave function becomes a real 
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exponential in the region R < Fa 2 /(a 2 — 1). As can be seen from (124|) . this is the region 
that is classically forbidden; the real, non-oscillatory, behaviour of the wave function is thus 
no surprise. 

Surprisingly, Eq. (I118p gives in fact already the complete class of solutions. This is shown 
in Appendix B. There thus exist no non-separating solutions, given the ansatz made where 
the full wave functional factorizes into functions on the respective lattice points. But we 
emphasize that we have succeeded in finding exact solutions to all quantum constraints. 
Other solutions to the full WDW equation and momentum constraints would necessarily 
couple the infinitely many shells comprising the dust cloud; to find them would demand a 
regularization scheme that is much beyond the scope of this paper. In a sense, the factor 
ordering chosen here leads to quantums states for which the WKB form is “exact”, cf. the 


analogous situation with the models discussed in 19] and 20 


For definiteness we will consider in the following only the positive sign in front of ar in 
(11181) . For a = 0, which corresponds to E — oo, the solutions are particularly simple. They 
read 


IFin ( r , T, R, a = 0) = const. ± (R + Fln(F — R)) (121) 

ll'out( r , T, R, a — 0) = const. ± (R + Fln(R — F)) , (122) 


where ‘in’ refers to the region inside the horizon and ‘out’ to the region outside the horizon. 
For a = 1, which corresponds to the marginal model E = 0, we obtain the solutions 


(r, T, R, a — 1) = const. + r ± 2 \[F 
ll'out( r , r, R, a — 1) = const. + r ± 2 Vf 


VR — VF tanh 1 
VR — VF tanh -1 


F 


These solutions were already found in [9] and discussed further in [H . 
The full solutions for the marginal model read 




ho„t( T '- R ' r 1 = *n,/out[r] exp 


drr(r)W± 0 ut (r,fl,r,o=l) 


(123) 

(124) 


(125) 


The solutions dO and have to be matched at the horizon, that is, both the states 
and their //-derivatives should agree there. As can be seen from (11231) and (1124|) . however, 
the phases of the states diverge there. One thus has to perform an analytic continuation: 
We write R — F = eexp (itp), e > 0, and compare the states at (p = 7t/2, that is, at 
R = Rh = F + ie. One Ends from this comparison that the states are related as 


vI 'm[A Rh, T] = 


j/TinjT] 

*rout[r] 


7T 

ex P I =F- 


dr T(r) F(r) ) ,[r, R h , T] . (126) 


The states and To Ut can be set equal at the horizon if we exploit the freedom to choose 
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\&r[T]. One thus gets a relation between Tp in [r] and 'I'± out [r], 

*? ln [r] = exp jf° dr r(r)nO «foJr] (127) 

= exp (t^(F 2 (oo) - J+0))j + ou ,[r] (128) 

= exp(^M 2 ) ^ out [r], (129) 

On the lattice this reads 

^rin( r i> Fi) = exp (+vruy Ft) ^ out {Ti, Fi) . (130) 


One can easily check that with this choice the derivatives of the states at the horizon coincide 
as well. We also remark that the alternative choice R = F — ie would lead to a switch of 
sign in the exponent of (1129(1 . 

We mention that the factor acquired when crossing the horizon might be connected to the 


entropy of the black hole 19]. In fact, we recognize that the absolute value of the exponent 


is one quarter of the Bekenstein-Hawking entropy. 

The solutions (II 18(1 are all solutions that can occur on the lattice for the factor ordering 
(11171) . Here we want to extend these solutions to the continuum. The solutions on the lattice 
contain two free parameters, a* and bi , 


_ ibiaTi/2 A^(aiTi+f dR 


T, = e 


\/l a i^ : i 


(131) 


(We have for simplicity considered only one sign in the exponent.) We know that a* is 
connected with the local energy E via 2Ej + 1 = 1 /of. In general we have E = E(r), and 
thus it would be natural to demand that a = a(r). But with this explict dependence on r 
the momentum constraint would not be satisfied. The only possible way out of this is to use 
an implicit dependence a(F(r )) in order to get solutions fulfilling the momentum constraint. 
The expression 


Vk[r,i?,T] =e mf^F{r))V 


xexpji 


dr T 


a(F(r))r + 



a(F{r))F 

~F 


(132) 


is still a solution, since the Hamiltonian constraint does not contain a derivative with respect 
to T or F. Hence we arrive at a family of solutions containing two arbitrary functions a(x) 
and b(x). It is clear that a(x) is connected with E via 2 E(r) + 1 = l/a 2 (F(r)). One might 
in principle wish to construct wave packets by superposing wave functions with different a, 
that is, with different energies. However, since the factor ordering depends on E, this does 
not seem feasible. 


IV. DISCUSSION 

Let us first summarize the main results of our paper. As for the classical part, we have 
extended the canonical formalism to the non-marginal case. We have shown in particular 
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that (1601) holds for all values of E. We have succeeded in expressing the observables F, E , 
and To, which determine the classical model completely, in terms of the canonical variables. 
We have presented a new and improved method to cope with the boundary action in the 
canonical formalism. This renders the formalism more transparent. We have presented the 
general relation between dust proper time and Killing time and noted the similarity of dust 
time with the generalized Painleve-Gullstrand time discussed in 17 . 

As for the quantum part, we have presented a lattice regularization that correctly imple¬ 
ments the momentum constraint in the continuum limit. For a particular factor ordering 
(that, however, depends on E) we have succeeded to give exact solutions to all constraints. 
We have shown that, given the general ansatz for the wave functional on the lattice, these 
are the only solutions. This means that these are the only solutions for which the states 
describing the dust cloud factorize into infinitely many states corresponding to the various 
shells forming the cloud. We have also extended these solutions into the continuum. 

The discussion of the quantum LTB model is far from being exhausted. We thus want 
to conclude with an outlook on future work. First, we have mentioned in Sec. III.C a 
possible connection between the exponential factor in (11291) and the Bekenstein-Hawking 
entropy. Perhaps it will be possible to recover this entropy as an entanglement entropy from 
these quantum states. This could provide a first step towards its general understanding. 
Second, it would be of great importance to study the possible singularity avoidance of the 
quantum LTB model. This would include also an understanding of the role of the naked 
singularities in the classical model. We recall that singularity avoidance was a main feature 
of the quantization of dust shells |>!]. Singularity avoidance was also shown in various models 
of loop quantum gravity |d]. A third issue would thus be to develop the loop quantization of 
the LTB model and compare its results with the results_obtained from the WDW equation. 


Extending the discussion of Hawking radiation in llj, |l2j], we plan to derive the corre¬ 


sponding Bogoliubov coefficients from the exact quantum states found in this paper. This 
should yield in the appropriate limits the thermal Hawking spectrum plus quantum gravi¬ 
tational correction terms. This is currently under investigation. 
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APPENDIX A: SIMPLIFICATION OF THE HAMILTONIAN CONSTRAINT 


Here we show how it is possible to eliminate Pp from (|58l) by using (|59l) and to obtain 
the new, relatively simple, version (l60j) of the Hamiltonian constraint. 

The momentum constraint (159|) can be solved for Pr, which gives 


Pr = 



R' 


(Al) 
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Inserting this expression into the Hamiltonian constraint (T58|) and noting that L 2 = R' 2 IR — 
4JFP/, 2 gives 


‘2P t \J 7 ~ 


2 + r' 2 = 


TR' ARP’ 


R 


+ 


P' 


yp T 


TP') 


TL 2 + ARt'P^P t 


R' 


Squaring the resulting equation yields 


(A2) 

(A3) 


ar ,2 p 2 l 2 + ar ,2 p 2 t' 2 = r 2 l 4 + 8Rl 2 Tt'p(,p t + i6R 2 t ,2 p' 2 p 2 


(A4) 


At , 2 P 2 (R'-PL 2 ) 


The second term on the left cancels with the corresponding term on the right. After dividing 
by L 2 one has 


4P /2 P 2 + ARP 2 t’ 2 = T 2 (R' 2 /R - 4PPp 2 ) + 8RYt'P' v P t . (A5) 

v -V-' 

L 2 

If we combine the second term on the left-hand side with the second and third term on the 
right-hand side, it is possible to use the momentum constraint to eliminate Pp, 

r 2 /? /2 

4P ,2 P 2 = —-- AR (r'P T - TPp) 2 . (A6) 

R " -v-' 

Rl2P R 

This then yields the desired expression for (j60j) : 

Pt + FP 2 R - A=0. (AT) 


APPENDIX B: UNIQUENESS OF QUANTUM SOLUTIONS 


Here we shall demonstrate that the separating solutions found from (110711 and (110811 are 
unique. 

Equation (110711 may be solved by the ansatz 


= a R w = s -^>. 

The function rj has to fulfill the integrability condition 

sin rj cost] cos?/ 

- ^d R T = —d T V , 

which leads to 


d T r\ = — \/~R tan rjdnrj — OrVr . 


(Bl) 


(B2) 


(B3) 


21 
























Inserting ansatz (IBID into (11081) gives another equation for 77 , 

1 A 

- i=d r 7j + cot rjdnTj - d R In T + — = 0 . 

v J- 

So the ansatz A = R8r yields 

d T r] = VI cot + VI 3 r In /r . 

The integrability conditions (IB3D and (1B5D have to be consistent. This gives 
— VI tan tiOrT] — dnVI = VI cot 7]8r7] + VI3r In /i . 

By an elementary manipulation we obtain 

Or In (Vlfi tanr/j = 0 , 

and thus 


a(r ) 

tan 77 = —■=, sm 7 ] = 


a(r) 


— , COS 7] = - ; - . 

[iVI V^TJH V 0/2 + If 

These equations can be used to give expressions for d R r) and d T r], 

dRri= j^ dR (A\ 


nVI 


a 2 A IT \nT) 

fill r, 

d T r] = —— —d T a . 

a- + i-iJ- 

Reinserting (IB9D and (IB10D into (1B3[) or (IB5D leads to 

d T a = %d R ii - . 

^ 2 


(B4) 

(B5) 

(B 6 ) 

(B7) 

(B 8 ) 

(B9) 

(BIO) 

(Bll) 


Since a is a function only of r, and /! and T are functions only of R , the above equation 
requires a = const. Hence we have 

^d Rf i = \d R T , (B12) 


which yields 


h' J 


h = 


P 


Vl - a 2 T ’ 

where P and a = P/a are constants. Then the unique solutions of (11071) and (11081) are 


(B13) 


W = ar 


d R 


Vl - a 2 T 

I 


(B14) 


which is just (11171) . 
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